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ABSTRACT 

Gamma ray bursts (GRBs) have recently attracted much attention as a possible way 
to extend the Hubble diagram to very high redshift. However, the large scatter in 
their intrinsic properties prevents directly using them as distance indicator so that 
the hunt is open for a relation involving an observable property to standardize GRBs 
in the same way as the Phillips law makes it possible to use Type la Supernovae 
(SNela) as standardizable candles. We use here the data on the X-ray decay curve 
and spectral index of a sample of GRBs observed with the Swift satellite. These data 
are used as input to a Bayesian statistical analysis looking for a correlation between 
the X-ray luminosity Lx{Ta) and the time constant Ta of the afterglow curve. We 
find a linear relation between log [Lx{Ta)] and log [Ta/{1 + z)] with an intrinsic scatter 
o'int = 0.33 comparable to previously reported relations. Remarkably, both the slope 
and the intrinsic scatter are almost independent on the matter density ^Im and the 
constant equation of state w of the dark energy component thus suggesting that the 
circularity problem is alleviated for the Lx — Ta relation. 
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1 INTRODUCTION 

The high fluence values (from 10^^ to 10~^ erg/cm'^) and 
the enormous isotropic energy emitted (~ 10^'' — lO^^erg) at 
the peak in a single short pulse make Gamma Ray Bursts 
(hereafter GRBs) the most violent and energetic astrophys- 
ical phenomena. Notwithstanding the variety of their dif- 
ferent peculiarities, some common features may be identi- 
fied looking at their light curves. Although GRBs have been 
traditionally classified as short and long depending on Tgo 
being smaller or larger than 2 s (with Tgo the time over 
which from 5% to 95% of the prompt emission is released), 
a recent analysis by Donaghy et al. ((2006)) has shown that 
this criterion has to be revised. Indeed, the existence of an 
intermediated class of GRBs have also been studied (Nor- 
ris & BonneU (2006); Bernardini et al. (2007)). As a result, 
the long GRBs are now further classified as normal and low 
luminosity with the latter ones probably associated with Su- 
pernovae (Plan et al. (2006); Dainotti et al. (2007)). 

Notwitstanding this classification, two phases are 
clearly visible in the GRB lightcurve, namely the prompt 
emission, where most of the energy is released in the 7 - rays 



in only tens of seconds, and an afterglow lasting many hours 
after the initial bursts. Early observations in the X-rays 
typically started several hours after the prompt emission so 
that only the late phase of the light curve could be charac- 
terized. It was then found that a phenomenological power - 
law, f{t,u) oc f'^u"'^ with {a, (3) ~ (-1.4,0.9), provided 
a reasonable fit to the observed data (Piro 2001). However, 
the launch of the Swift satellite, whose aim is also to observe 
GRBs X-ray (0.2-lOkeV) and optical (1700-6500 A) after- 
glows starting few seconds after the trigger, revealed a more 
complex behaviour. The soft X - ray light curves must indeed 
be divided in two different classes (Chincarini et al. 2005) 
according to the steep or mild initial decay. Most of the 
observed GRB afterglows belong to the first group, show- 
ing what has been called a canonical behavior (Nousek et 
al. 2006) described by a broken power - law. After the ini- 
tial steep decay (with slope 3 < Qi < 5), the light curve 
shows a shallow decay (0.5 < Q2 < 1) followed by a some- 
what steeper decay (1 < Q3 < 1.5) beyond 2x10* s. These 
power - law segments are separated by two corresponding 
break times with tbi < 500 s and 10^ s < tt,2 < lO" s. A 
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new systematic study using GRBs observed with XRT re- 
veals a still more complex behavior with diflFerent power - law 
slopes and break times (O' Brien et al. 2006; Sakamoto et 
al. 2007). A significant step forward has been represented by 
the analysis of the X - ray afterglow curves of the full sample 
of Swift GRBs showing that all of them may be fitted by the 
same analytical expression (Willingale et al. 2007). 

Finding out a universal feature for GRBs is the first im- 
portant step towards their use as distance indicator. To this 
aim, one has indeed to look for a universal relation linking 
observable GRBs properties so that their intrinsic luminos- 
ity may be estimated from directly measurable quantities. 
Previous attempts along this road arc represented by the 
EiSO — Epeak (Amati et al. 2002), E.^ — Epeak (Ghirlanda 
et al. 2004; Ghirlanda et al. 2006), L-Ep^ak (Schaefer 2003), 
L — Tiag (Norris et al. 2000), L — V (Fenimore & Ramirez - 
Ruiz 2000; Riechart et al. 2001), L — trt- Moreover, three - 
parameter relations have also been proposed such as, e.g., 
the EiSO — Ep — tt (Liang & Zliang(2005)) and that proposed 
by Firmani and collaborators (Firmani et al. 2005, 2006). 
On the other hand, some attempts have also been made to 
compare these empirical correlations with the model depen- 
dent ones (Nava et al. (2006); Guida et al. 2008). The above 
quoted two parameters correlations have then been used by 
Schaefer (2007, hereafter S07) to construct the first reliable 
GRBs Hubble diagram extending up to a ~ 6 opening the 
way towards the use of GRBs as cosmological probes (see, 
e.g., CapozzicUo & Izzo 2008 and rofs. therein). 

In this letter, we present a possible alternative route 
towards standardizing GRBs as distance indicator. To this 
aim, we use the data in Willingale et al. (2007) to look for 
a possible correlation between the X-ray luminosity at the 
break time Ta and the Ta itself. The data used are presented 
in Sect. 2, while Sect. 3 deals with the statistical tools and 
the results. Conclusions are summarized in Sect. 4. 



2 THE DATA 

Willingale et al. (2007, hereafter W07) have examined the 
X - ray decay curves of all the GRBs measured by the Swift 
satellite then available. Their analysis shows that all of them 
may be well fitted by a simple two components formula, 

namely : 



f{t)= f pit) + fait) 



(1) 



where the first term accounts for the prompt 7 - ray emission 
and the initial X - ray decay, while the second one describes 
the afterglow. Both components are given by the same func- 
tional expression : 

Fc exp (^ac - exp (^y) for t <Tc 

hit) = { ' (2) 

Fe(^)""%xp(-|) for t>T. 



where the transition from the exponential to the power, law 
decay takes place at the point (Tc, where the two func- 
tional sections match in value and gradient. The parameter 
Oc determines both the time constant of the exponential de- 
cay (given by Tc/oc) and the slope of the following decay, 
while tc marks th e initial r ise and the time of maximum flux 
occuring at t = ^JtcTc|a.c■ Denoting with the suffix p and o 



quantities for the prompt and afterglow components, Eq.(2) 
may be inserted into Eq.(l) to give an eight parameters ex- 
pression that can be fitted to the X - ray decay curve in or- 
der to both validate this expression and determine, for each 
GRB, the corresponding parameters. Such a task has been 
indeed performed by W07 using all the 107 GRBs detected 
by both BAT and the XRT on Swift up to August 1st 2006. 
The fit procedure and the detailed analysis of the results are 
presented in W07, while here we only remind that the usual 
fitting in the log (flux) vs log (time) provide estimates 
and uncertainties on the time parameters (log Tp, log Ta) and 
the products (logFpTp, logF^Ta). 

W07 also performed spectral fitting with XSPEC (Ar- 
naud 1996) to BAT (for the prompt phase) and XRT (for 
later phases) data to estimate the spectral index during dif- 
ferent phases. Due to the limited frequency range, the GRB 
spectrum may be simply described by a single power - law, 
$(-E) oc E'^ , with the slope (3 depending on the time when 
the spectrum is observed. W07 reported four different values 
of namely [3p (for the prompt phase), I3pd for the prompt 
decay, /3o for the plateau observed at the time Ta, and /3ad 
for the afterglow at t > Ta. Actually, the data coverage is 
not sufficient to measure all of them for the full sample so 
that, for the weakest bursts, only Pp and Ppd are available. 
Provided /3 is known, it is possible to estimate the GRB 
luminosity at a given time t as : 



Exit) =4TvDl{z)Fxit) 



(3) 



where DLiz) is the luminosity distance at the GRB redshift 
z, and Exit) is the fiux (in erg/cm^ /s) at the time t, K- 
corrected (Bloom et al. 2001) as : 



Exit) = fit) X 



j^-- E^E)dE 



(4) 



with (E,n„-i, E,nax) ~ (0.3,10) kcV set by the instrument 
bandpass. Note that Eq.(3) is the same as Eq.(8) in S07 the 
only difference being the integration limits of the integral 
at the numerator. Actually, while S07 is interested to the 
bolometric luminosity, we are here concerned with the X- 
ray one so that we integrate only over this energy range. 

Using the data in W07, we compute the X-ray lumi- 
nosity at the time Ta so that we have to set fit) = /(Ta) 
and f} — fSa in Eqs.(3) and (4). Actually, rather than us- 
ing Eq.(l), wc set /(T^) = faiTa) since the contribution of 
the prompt component is typically smaller than 5%, much 
lower than the statistical uncertainty on faiTa). Neglecting 
/p(T'a) thus allows to reduce the error on FxiTa) without 
introducing any bias. This latter error is then estimated by 
simply propagating those on [3a, logTa and logFaT^ thus 
implicitly assuming that their covariance is null*. Should 
this not be the case, we are underestimating the final error 
on Lx(Ta). We have, however, checked that our main results 
are unaffected by a reasonable increase of the errors. 



* Note that the covariance matrix is not reported in W07, where 
the parameters of interest arc given with tlieir 90% confidence 
ranges. Following Willingale (priv. comm.), we have assumed 
independent Gaussian errors and obtained la uncertainties by 
roughly dividing by 1.65 the 90% errors. Moreover, we prelimi- 
nary correct for asymmetric errors on logFoTo and log Ta (when 
present) following the prescriptions in D' Agostini (2004). 
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As a final important remark, we note that the presence 
of the luminosity distance Dl{z) = {c/Ho)dL{z) in Eq.(3) 
constrains us to adopt a cosmological model to compute 
LxiTa). We use a flat ACDM model so that the Hubble 
free luminosity distance reads : 



di(«) = (l + ; 



dz' 



^/f^M{l + z'y' + {l-nM) 



(5) 



In agreement with the WMAP five year results (Dunkley et 
al. 2008), we set (Qm, ft) = (0.291, 0.697) with h the Hubble 
constant Ho in units of 100 km/s/Mpc. 
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3 A LUMINOSITY -TIME CORRELATION 

In order to standardize GRBs to use them as possible dis- 
tance indicator, we need to find a correlation between the 
luminosity and a directly observable quantity. Should such a 
relation be found, one can then use the observed fiux and the 
estimated Lx to infer Dl{z) and then construct the GRB 
Hubble diagram. Let us suppose that a power - law relation 
exists between two quantities R and Q as R = AQ^. In log- 
arithmic units, this roads log R — a + b log Q with a — log A 
and b = B. Typically, both R and Q will be known with mea- 
surement errors (an, <tq) so that the statistical uncertainties 
on (logJi,logQ) will be given by (cth/T?, crg/Q) x (1/ In 10) 
respectively. These errors may be comparable so that it is 
not possible to decide what is the independent variable to be 
used in the usual fitting analysis. Moreover, the relation 
R = AQ^ may be affected by an intrinsic scatter aint of un- 
known nature that has to be taken into account. In order to 
determine the parameters (a, b, aint), we can then follow a 
Bayesian approach (D' Agostini 2005) thus maximizing the 
likelihood function C{a, b, aint) = exp [—L{a, b, aint)] with : 



L{a,b,ai„t) 



2 I 2 I 7 2 2 
int + '^■y,- + 



lE'nia 

1 \ ^ (vi-a- bxj)'^ 

2 2^ al, +al+b'^al 



(6) 



with {xi,yi) = (log Qi, log _Ri) and the sum is over the J\f 
objects in the sample. Note that, actually, this maximization 
is performed in the two parameter space (6, aint) since a may 
be estimated analytically as : 



t/i - bxi 



+ al + V^c 



+ al. + b'^c 



(7) 



so that we will not consider it anymore as a fit parameter. 

We use this general recipe to look for a correlation 
between the X-ray luminosity (in erg s~^) at the time 
Ta and Ta (in s) itself, i.e. we set y = log [Lx (To)] and 
X = log [Ta/ (1 + z)], where we divide time by (1 -|- z) to ac- 
count for the cosmological time dilation. Note that, since f3a 
is needed to compute Lx(Ta), we have to reject most of the 
107 GRBs reported in W07 because this is not known. We 
thus end up with a sample contanining A/" = 32 with both 
log [Lx{Ta)] and log [Ta/{1 + z)] measured^. The Spearman 
rank correlation turns out to be r = —0.74 suggesting that a 

t ASCII tables with all the quantities needed for the analysis and 
the Mathematica codes used are available on request. 



log T./(l+z) 



Figure 1. Best fit curves superimposed to the data with the 
solid and dashed lines referring to the results obtained with the 
Bayesian and Levemberg - Marquardt estimator respectively. 



power- law relation between Lx{Ta) and Ta/(1 + z) indeed 
exists thus motivating further analysis. 

We then apply the maximum likelihood estimator de- 
scribed above in order to determine both the slope and the 
intrinsic scatter of the Lx — Ta correlation thus finding out : 

(o, 6, aint) = (48.54, -0.74, 0.43) . 

Defining the best fit residuals as 5 = yobs — Vfit, we can 
qualitatively estimate the goodness of the fit by consider- 
ing the median and root mean square which turn out to be 
{S) = —0.08 and Srms = 0.52 indeed quite small if compared 
to the typical log [Lx (Ta)] values. It is also worth noting that 
5 does not correlate with the other parameters of the fit flux, 
while the value r = —0.23 between S and z favours no signi- 
ficative evolution of the Lx — Ta relation with the redshift. 
The best fit relation is superimposed to the data in Fig. 1 
where we also present the best fit obtained by the usual 
fitting technique. In this case, the best fit parameters are 
obtained by minimizing (through a Levemberg - Marquardt 
algorithm with 1.5cr outliers rejection) a merit function 
given by the second term in Eq.(6) with ay^ = aint ~ 0, i.e. 
we (erroneously) assume that there is no scatter and that 
the errors on \og[Lx{Ta)] are negligible. This alternative 
method gives as best fit parameters : 

(a.b) = (48.58, -0.79) 

in good agreement with the above mgiximum likelihood esti- 
mator so that we argue that our results are independent on 
the fitting method. However, since the Bayesian approach is 
better motivated and also allows for an intrinsic scatter, we 
hereafter eligc this as our preferred technique. 

In an attempt to reduce the intrinsic scatter in the 
above correlation, we have analysed the best fit residu- 
als noting that the higher ones are obtained for GRBs 
with luminosities smaller than lO**^ erg and time param- 
eter log [To/(l + z)] > 5. We therefore repeated the above 

analysis using only 28 out of 32 GRBst satisfying the two se- 



t The four GRBs excluded are: GRB050824, GRB060115, 
GRB060607A and GRB060614. While the first two appear to be 
unaffected by any problem, for the latter two, the data cover less 
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log T,/(l+z) 

Figure 2. Same as Fig. 3, but only using GRBs with 1 < 
log [Ta/{1 + z)\<b and log [Lx{Ta)\ > 45. 



lection criteria 1 < log [Ta/(1 + ^)] < 5 and [Lx(Ta)] > 45. 
Using the maximum likelihood estimator, we get : 

(a, b, aint) = (48.09, -0.58, 0.33) 

with {6} = —0.06 and 6rms = 0.43. The reduced intrinsic 
scatter and the smaller fit residuals suggest us that, what- 
ever is the unknown mechanism originating the Lx — Ta 
relation, this is better effective for the class of GRBs sat- 
isfying the above selection criteria. The data and the best 
fit curve are shown in Fig. 2 where the dashed line refers 
to the results obtained with the minimization giving 
(a, b) = (48.07, —0.60) reported here for completeness. 

The Bayesian approach used here also allows us to quan- 
tify the uncertainties on the fit parameters. To this aim, for 
a given parameter pi , wc first compute the marginalized like- 
lihood Ci{pi) by integrating over the other parameter. The 
median value for the parameter pi is then found by solving : 



rv 



.,med ^ 

Ci{pi)dpi = - 



Ci{pi)dpi 



(8) 



The 68% (95%) confidence range {pi,i,Pi,h) are then found 
by solving : 



Ci{pi)dpi 



£.i{pi)dpi = — — 



Ci{jpi)dpi 



Ci(pi)dpi 



(9) 



(10) 



with e = 0.68 (0.95) for the 68% (95%) range respectively. 

For the fit to the full datasct, we get : 



r,.+0.20 -HO. 41 

-U. ('t.Q.ig -0.39 



_ a /IS+015 -HO. 35 
CTint — 0.'±0_Q iQ _o.l8 J 



while it is : 



h- -n t:«+018 +0-38 
O — U.00_o.i8 -0.37 



_ n qQ-l-014 -HO. 33 
CTint — 0.0»_Q 11 -0.20 



for the selected subsample. 



4 DISCUSSION AND CONCLUSION 

The high Spearman correlation coefficient, the low value of 
the fit residuals and the modest intrinsic scatter renders the 
Lx — Ta relation presented above a new valid tool to stan- 
dardize GRBs. It is worth stressing that Lx — Ta needs only 
two parameters and one of them is directly inferred form 
the observations minimizing the effects of the systcmatics 
errors. Furthermore the redshift range covered is large ex- 
tending from 0.54 (0.125) up to 6.6 for the selected (full) 
sample far beyond the maximum redshift affordable with 
Type la SNc [z ~ 1.7). Should this correlation bo confirmed 
by future higher quality data, one could then combine it with 
the other relations yet available in literature to work out a 
GRBs Hubble diagram deep into the matter dominated era 
thus representing an outstanding cosmological test. 

To this end, it is worth comparing the Lx — Ta relation 
with other ones quoted in literature. When performing such 
a comparison, however, one should take into account the 
differences in the cosmological model adopted and the fitting 
method used. In particular, the choice of how the best fit 
parameters are estimated may have an important impact on 
the estimate of the intrinsic scatter with the usual fitting 
leading to an underestimate of (Jint- On the other hand, 
changing Q,m in the framework of the fiat ACDM scenario 
have a profound impact on Oint with higher giving rise 
to lower (Tint values (Basilakos & Perivolaxopoulos 2008). In 
order to account for both these issues, one should therefore 
test all the above correlations using the same statistical tools 
and cosmological model, a task wc will address elsewhere. 

As is well known, the paucity of local (i.e., z < 0.1) 
GRBs represents a serious problem for any attempt to stan- 
dardize GRBs since it is very difficult to directly calibrate 
any relation. This problem may be partly overcome by fit- 
ting the correlation in a subsample of GRBs lying at similar 
redshift (Amati 2008). However, as a general rule, in order 
to evaluate the GRB luminosity, a cosmological model has 
to be adopted thus leading to the circularity problem. Al- 
though addressing this problem in detail will be the subject 
of a forthcoming work, wo have here investigated what is 
the effect of changing the cosmological model by using our 
maximum likelihood estimator to determine the parameters 
{a,b,Gint) as function of Q.m in a flat ACDM modcl^. Wc 
find the remarkable result that both the best fit parame- 
ters (6, (Tint) and the rms residual 5rma are almost insensi- 
tive to the value of Q^m- Indeed, b runs from b ~ —0.590 
to 6 ~ —0.565, while aint increases from (Tint — 0.335 to 
aint — 0.340 for Q,m going from 0.2 to 1.0. As a further test, 
we generalize the ACDM model varying not only the matter 
density parameter Q.m, but also the equation of state w of 
the dark energy component (with w = —1 for the ACDM 
model). For —1.3 < w < —0.7, neither b nor rjint signif- 
icantly change confirming the qualitative results obtained 
for the ACDM scenario. Although a more detailed analysis 
is needed, we therefore argue that the circularity problem is 
alleviated by the use of our Lx — Ta relation. 

The encouraging results discussed above are serious ar- 
guments in favour of the Lx — Ta relation as a further tool 



than 50% of the Tgo range . Moreover, for GRB060607A, the 
prompt component dominates over the afterglow one so that our 
approximation f{Ta) ^ fa{Ta) is not valid anymore. 



§ To be precise, we let CIm running from 0.2 to 1 and adjust h 
so that Q.Mh'^ is fixed to the same value adopted above. 
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towards the standardization of GRBs as distance indicator. 
Should these first evidences be furtherly enforced by more 
data, the combined use of full set of GRBs correlations dis- 
covered insofar could opened the road towards making GRBs 
the high rcdshift analog of SNela as cosmological probes in 
the not too distant future. 
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